High-spin polaron in lightly doped Cu02 planes 
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We derive and investigate numerically a minimal yet detailed spin polaron model that describes 
lightly doped Cu02 layers. The low-energy physics of a hole is studied by total-spin-resolved exact 
diagonalization on clusters of up to 32 Cu02 unit cells, revealing features missed by previous studies. 
In particular, spin-polaron states with total spin 3/2 are the lowest eigenstates in some regions of 
the Brillouin zone. In these regions, and also at other points, the quasiparticle weight is identically 
zero indicating orthogonal states to those represented in the one electron Green's function. This 
highlights the importance of the proper treatment of spin fluctuations in the many-body background. 

PACS numbers: 71.10.Fd,75.10.Jm,71.38.-k,74.72.-h 



Introduction: A full understanding of the physics of 
a Cu02 layer doped with a few holes has still not been 
achieved, despite continuous effort . Recent high res- 
olution angular resolved photoemission (ARPES) studies 
0t5i| on the insulating charge-transfer gap parent com- 
pounds reveal major puzzles: do quasiparticles of one 
electron nature exist and if so what is their energy and 
momentum? Why are the first visible electron removal 
states so broad, of 300 meV at 300K and decreasing lin- 
early with temperature T, and what causes the very ap- 
parent T-dependent change of line shapes? Is the mo- 
mentum dependence of the lowest energy structure re- 
lated to the pseudo-gap formation at higher hole concen- 
trations? Recent neutron experiments performed in the 
pseudo-gap phase reported magnetic response through- 
out the Brillouin zone, not restricted to the region of the 
much discussed magnetic resonance Q • These and other 
issues including the broken local 4-fold symmetry, which 
is taken for granted in single-band models, seen in scan- 
ning tunneling probe (STM)0 and X-ray scattering Q 
remain either open questions or are controversial. 

It is widely believed that a complete description of a 
single hole in a spin-i 2D antiferromagnet (AFM) with 
full quantum fluctuations could provide the answers to 
these questions, as well as clues for understanding the 
origin of the non- Fermi-liquid behavior and the supercon- 
ducting ground-state in the higher hole density region. 
Of course, consideration of exotic many-body scenarios 
[lo| . or of coupling to lattice vibrations [Ilj are excit- 
ing developments; however, a detailed modeling of the 
hole-f AFM is a crucial first step to understand the sig- 
nificance of such additions. This problem is very difficult 
because of the complicated nature of the 2D AFM back- 
ground, whose quantum fluctuations in the presence of 
doped holes were never fully captured for a large Cu02 
lattice. Recent technical developments [l^ allow us to 
present the first such results for samples with up to 32 
Cu and 64 O, in this Letter. Our work also reveals impor- 
tant failings of the single-band models used extensi vely 
in the literature, and which we briefly review below [13|. 

Microscopic hole- AFM interactions have been studied 



in models with one U^UM^ t^o U^^, three |2l|, or 
more 22, [2^ bands. While exact analytical solutions 
seem to be out of reach, numerical studies are always 
carried out with compromises such as the use of small 
clusters and variational approaches [l^. Given these dif- 
ficulties and the drive to find the simplest model, the one- 
band models are unsurprisingly the most studied 114 1. 
While certain higher-energy aspects observed by XAS, 
EELS and STM 0, ^ cannot be described using one 
band, the significance of omitting other bands in the low- 
energy scale cannot be quantified without a comparison 
to unbiased solutions of more detailed models. 

Cuprates exhibit charge-transfer band-gap behavior 
with mobile holes located mainly on anion ligands and 
unpaired electrons on cation d-orbitals @. One-band 
models use superexchange [i^ and Zhang-Rice singlets 
(ZRS) [m to reduce the ( — n)-clcctron problem to one 
of n holes in an AFM background, often modeled as a 
Neel background with spin-waves. To reach agreement 
with experiments, such models must be tweaked at least 
by adding longer-range hopping [i^, 27 1. One trade-off 



for their elegance is the use of momentum-independent ef- 
fective parameters, even though it is well known that the 
ZRS state has a strong k-dependent renormalization [l^ . 
The impact of such approximations must be verified for 
all k with models that distinguish anion and cation sites. 

In this Letter we study a single hole in a model that 
includes the O 2p orbitals explicitly and full quantum 
fluctuations of the AFM background. This results in 
spin-polaron solutions absent from other models or ap- 
proximations used to date. The energy dispersion of the 
lowest energy electron removal state is similar to that of 
the ZRS; which effectively locks the O hole in a singlet 
with one of the two adjacent copper sites. Without such 
restrictions, our resulting wavefunction features the hole 
forming a stable S = ^ three-spin polaron (3SP) with 
its two neighbour copper sites. In some regions of the 
Brillouin zone, a 5 = 1 quantum fluctuation binds to the 
S = polaron yielding a low-energy 5' = | state invisible 
at T = to ARPES. For all momenta, the 5' = | states 
are found to be within ^ J/2 of the S = band. These 
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FIG. 1. (a) Two adjacent unit cells of the Cu02 plane. The 
orbitals kept in the 3-band model of Eq. ([T]) are shown, with 
white/shaded for positive/negative signs. The two e vectors 
(solid arrow) and the four 5 vectors (dashed arrow) are also 
shown, (b) Sketch of a virtual process of Tswap. 
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and other results inconsistent with one-band models, and 
their experimental implications, are discussed below. 

We start with the three-band p ~ d model which ex- 
hibits the basic physics of a hole doped charge transfer 



gap and insulating spin 1/2 antiferromagnet [21 



HsB — Tpd + Tpp + Apd ^ nij^^^cr 

+ Upp ^ ,|71;+e,4. + Udd ^ '^;,t"i,4- (1) 



where ni^cr = dj^^di^a, «i+e,o- = Pi+^^crPi+>^,<r count holes 
in Cu ^d^2_y2, respectively O 2p^/y orbitals (see Fig. 
1) and Udd > Upp > Apd describe Hubbard and charge 
transfer interactions. Nearest neighbor (NN) Cu-0 hop- 
ping Tpd = tpdJ2[ipl+e,a-pl-e,a)d'i:'^+f^-^-] '^^ included, as 



is hopping Tj 



pp — i-pp ^ •3di^i+e+5_cr-Pi+e,o- "pp . 

Pi+3e a)Pi+<^,'^ between NN and certain NNN O sites. For 
NN hopping by (5 = {5x,8y), ss = SxSy/\Sj:Sy\. 

In a half-filled, laige-U system with no hopping^ the 
ground-state (GS) has a hole at each Cu site: Yl ^1 ^JO) ~ 
YiWi)^ with the usual 2^ spin degeneracy. An electron 
removal adds a hole in an O orbital, so the doped GS 
is cr n I '^ith X 2^+^ degeneracy. We study 
the behavior of such anion holes when the hopping is 
turned on, in the framework of superexchange. The idea 
is reminiscent of studies such as Refs. 3, 1^ 22 1; how- 
ever, these also used further approximations. A detailed 
comparison of our Hamiltonian versus those used in these 
references is provided in the Supplementary Material [l3| . 

Model: Noting that all Tpd processes increase energy 
by either U and/or Apd, we derive the effective model for 
the states pj^^ aHl'^i) ^'^ fisj : 

Hcs = Tpp + Tswap + Hj^^ + Hj^^ (2) 

where the 0-0 hopping of the hole is supplemented by: 

Tswap = ~tsw^S^pl^^^^^^^pi+,^^,\(jlJ{(Ti^J (3) 

Hjpa = Jpd X! ' '^'±'= (^) 
Hj^ia = Jdd ^ Sl±2e ■ S'in^(l - ?^/±e,<T) (5) 

Using tpd = l.2>eV, tpp = 0.65eT^, t'pp = 0.58tpp, Apd = 
Z.&eV , and Upp = AeV 14|, we scale the parameters in 



FIG. 2. a) Energy and b) quasiparticle weight (bottom) for 
the lowest eigenstates with St = \ and | vs. momentum. 
Different sets are shifted so as to have the same GS energy. 



units of Jdd to find their dimensionless values to be tpp = 
4.13, tsu, = 2.98, and Jpd = 2.83. 

While we find that the 3-spin polaron (3SP) [l^ plays 
an important role, our approach is different from previ- 
ous work [l3| by recognizing (i) Tpp's role as a coherence 
facilitator rather than a mere correction; (ii) its comple- 
menting process Tswap, illustrated in Fig. 1(b), (iii) sup- 
pression of superexch ang e along the bond inhabited by 
the hole, see Eq. ^ [13[, and iv) total-spin (St) eigen- 
states are studied explicitly. We push the computational 
limit to perform S'r-resolved exact diagonalization (ED) 
of a topologically superior [i^l cluster of = 32 Cu02 
unit cells, treating the AFM background exactly. ED 
provides the transparency, flexibility, and neutrality to 
support new results. The price for a systematic mapping 
of the excited states is the limited k resolution. A^ = 16 
results are provided to check for size dependence. 

We find that all low-energy eigenstates have a total 
spin of either St = ^ or |. The z-projections for each 



St are degenerate. The St = ^ subspace is due to the 
s = i hole mixing with various = background states, 
including the AFM GS, or mixing with the S = 1 back- 
ground states, including the "single-magnon" states. The 
s = i carrier can also mix with = 1 or 2 background 
states to yield the St = § subspace; we explicitly con- 
sider such states here for the first time. The partition of 
the S^ subspace into separate St sectors was managed 
by the optimizations of Ref. [l2j . Unlike there, no basis 
truncation was employed here for rigorous results. The 
(k, St = ^, S^ = ^) sector contains ^0.44x10^ states. 

Results: Fig. [2ja) shows the lowest eigenenergies. The 
GS has k = (f , f ) and 5t = ^, is consistent with the 
3SP but can also be thought of in terms of ZRS [isj . 
Remarkably, we find similar dispersion along (0,0)^(7r,7r) 
and (0,7r)— >(7r,0) without having to add longer range hop- 
ping or fine-tune parameters as is needed in one band 
models. The biggest surprise, though, are the low-lying 
St = ^ states which go below the ^ states near (0,0) 
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FIG. 3. {Cx{S,a)) for the lowest energy state at (a) {^,^) 
with St = |, and (b) at (7r,7r) with St = |. The darkly- 
shaded bullet denotes the oxygen position at I + Cx- Each 
bullet shows the correlation value between the two sandwich- 
ing Cu sites. The central 12 Cu sites are shown; the corre- 
lations between the other 20 Cu spins converge fast towards 
the AFM value of ~-0.33. {Cy{5,a)) is the P^^y reflection. 



and (tt.tt). Finite-size analysis, discussed in the Supple- 
mentary Material, reveals that the S't = | states are 
stable polarons at least in the regions marked by thick 
solid lines in Fig. ^a). Thus, a S't = 5 quasiparti- 
cle cannot describe the low-energy states throughout the 
BZ. To compare the lowest energy states on both sides of 
the crossing, we note that the lowest = ky eigenstates 
have odd parity upon a Px^j, reflection (FigU^) so they 

can be expressed as + Px^y)J2 '^''''''pl+e^,aWJ)x- 

The band-crossing results in noticeable change in the ex- 
pectation values of the correlation function: 



CxiS, a) =2 ^ Si+s ■ Si+s+ 



(6) 



which measures the correlation between two neighboring 
Cu spins at a distance S from the hole. (C) ranges from 
-3/4 for singlet, to —0.33 for 2D AFM, to 1/4 for triplet. 

Fig. 13^ shows (Cx) when the hole is located at the 
darkly shaded bullet, in the GS: k = (f,f), St 
{(Cy) is a reflection with Px^^y for kx 
affects the AFM order in its vicinity. Because of the 
hole-spin exchange Hj^^ and the blocked superexchange 
between the two Cu spins neighboring the hole, these 
"central" spins have triplet correlations, of ^0.13. Also, 
(TJjpj) ^ -0.9Jpd, showing that locally this is consistent 
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ky). The hole 



with the 3SP solution [13| . More interesting are the cor- 
relations with the other 3 neighbors of each of these cen- 
tral Cu spins: with two of them, there are robust AFM 
correlations of ~ —0.22, while with the third the cor- 
relation nearly vanishes (lightly shaded bullet). This is 
counterintuitive if one views the system as a fluctuat- 
ing Neel background, where a spin-flip would change the 
spin-spin correlation to all four neighbors. Although the 
two central Cu spins have | weight in triplet configu- 
ration which is hardly bi-partite, long-range AFM order 
cannot be automatically discounted |29| . Indeed, the cor- 
relations we find are consistent with such order, except 
for the zigzag of 3 bonds shown by shaded bullets. This 
strange shape is dictated by the hopping mechanism. For 
a Bloch wave, 0-0 hole hopping in the upper- left/lower- 
right direction yields a phase shift of e'°/e'(*''==~''!'^ and 





FIG. 4. {Cxjy{5, a)) for lowest state at k=(0,7r) with 5t = 1/2. 



hence constructive interference if kx = ky. In contrast, 
hopping in the upper-right/lower-left direction yields a 
phase shift of e*'==/e and the interference is scaled 
down by cos{kx/y). In the GS, having a mixture of sin- 
glets and triplets upper-left/lower-right to the O hole 
lowers energy with the least disturbance to AFM order. 
Thus, the two outside zigzag bonds are triplet "distur- 
bance tails" pointing orthogonal to the momentum direc- 
tion. This is very different from the ZRS, which freezes a 
Cu spin by intraplaquette coherence with its four O sites. 

Fig. shows the correlation values when the S't = § 
polaron becomes the lowest energy state at {it, it). The 
results look similar at (0,0). (Hj^^) remains ^ — 0.9 Jpd, 
but there are now four more heavily disturbed bonds. 
This further supports this being a stable polaron with 
an extra magnon bound locally close to the O hole. We 
stress here that this | polaron is formed by a spin dis- 
turbance around the 3SP. This is very different from the 
S = I excitation local to Hj^^ with energy p^ . 

Fig. [2b shows the quasiparticle weight Z(k) for the first 
electron removal state. The major difference from other 
models is that 2'(k) = in three regions: a) Z(0,0) = 
Z{tt,7t)=Q because here the lowest eigenstate has St = 
I which due to spin-conservation is not in the Krylov 
space of any St = 5 state [13], and b) Z(0,7r)=0 even 
though this is a St = 5 state (see below). The t-t'- 
t"-J model treatment does not conserve St, resulting in 
Z(0, 0) ~ 0.1 and a finite ^(tt.tt) Our Z{k) is smaller 
everywhere than that of the t-t'-t"-J model, suggesting 
less "free particle" nature of the polaron. 

The lowest energy state at k = (0, tt) has St = 5, but 
its Z = because the state is not in the Krylov space 
of an electron removal [l^. This seems to be due to 
symmetry, although we do not yet fully understand this. 
States with finite but small Z are at least 0.006J higher 
in energy. Their close existence above the Z = lowest 
state may be related to pseudogap phenomena in this 
region; however, this needs to be investigated in more 
detail. Fig. |4] shows the correlation for this Z = lowest 
state. Compared to the GS (see Fig. [3^), there are 2 
more disturbed bonds as required by the reflection parity 
about k. This larger disturbance range is accompanied 
by more negative (AFM) correlation values. 

Although we are restricted to rather low momentum 

resolution, more can be said about the Es — Ei = 

2 , 2 

band crossings in the nodal direction by looking at the 
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k points between which the difference switch signs. The 
observation in Fig. [2^ is that, going away from the (f , f) 
GS, — El is larger towards (0, 0) than towards (tt, tt). 
The band crossing would induce an abrupt change 

in Z(k) from non-zero to exactly zero, irrespective of the 
Z[k) value on the finite side. The larger E^ towards 
k ~ (0, 0) suggests that the non-zero region extends more 
towards k = (0,0) than towards k = {tt,tt). Fig. [5^ also 
shows that the | states get pushed further down a.s N ^ 
oo so the crossing is expected to be closer to the ^) GS. 
This is consistent with ARPES which indeed observed an 
abrupt peak suppression in the nodal direction as well as 
the peaks surviving longer towards fc = (0, 0) 

Even when the S't = | states are not lowest in energy, 
they hug the S't = ^ band. This provides a ^ Jdd/2 en- 
ergy scale for spin excitations. At finite T, as magnons 
become thermally activated, these | states become "visi- 
ble" to ARPES. This suggests a T-dcpendent broadening 
mechanism of < scale, which, coincidentally, is the 

same energy scale recently linked to phonons, as another 



in the anti-nodal region is still being investigated. 
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possible source for this broadening 

Recent neutron experiments on samples at higher dop- 
ing reveal ^ 50meV magnetic response centered at q=0, 
away the AFM resonance momentum 01 . The bottom of 
the single-particle band structure in Fig 2a) indeed has 
a q=0 l/2-to-3/2 excitation of this energy scale. Our re- 
sults so far are restricted to a single hole; nevertheless, it 
has been pointed out that the q=0 magnetic excitation 
can be explained by involving spins on oxygen sites [30| . 
which certainly are present in our results. 

In addition to the low-energy 3/2 polaron band, there 
are internal energy scales of the local 3SP since Hjp^ 
also has a S" = ^ doublet and a 5 = | quartet separated 
in energy by Jpj^ and 3 Jp^/S from the lowest energy state. 
Magnetic excitations of these energy scales have been ob- 
served via inelastic resonant X ray scattering even for 
doped samples without long-range AFM order (sij . 

Summary: We solved a detailed model which includes 
the O sites and takes full account of the AFM quantum 
fluctuations, for large = 32 clusters. The phases of 
the p and d orbitals lead to phase coherence via Tpp -|- 
Tswap This is re-enforced by Hj^^ and the blocking 
of the AFM superexchange, making corrections such as 
Txondo negligible. While the dispersion is similar to that 
measured by ARPES without any fine-tuning, the lifting 
of Cu-0 singlet restriction present in ZRS-based models 
leads to wavefunctions of different nature, i.e. the 3SP 
where the O hole correlates with both its neighbour Cu 
sites. This model also provides low-energy channels for 
S = 1 excitations. Z{k) was found to be identically zero 
in certain regions of the BZ for two reasons: 1) the spin- 
I of the lowest energy state close to (0,0) and (7r,7r); 
and 2) around the antinodal region because of the lowest 
energy state there being exactly orthogonal to the single 
electron removal state. The detailed nature of the state 
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